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We report on new observations of the motion of zodiacal 
dust using optical absorption line spectroscopy of zodi- 
acal light. We have measured the change in the profile 
shape of the scattered solar Mg I A5 1 84 line toward sev- 
eral lines of sight in the ecliptic plane as well as the eclip- 
tic pole. The variation in line centroid and line width as a 
function of helio-ecliptic longitude show a clear prograde 
signature and suggest that significant fraction of the dust 
follows non-circular orbits that are not confined to the 
ecliptic plane. When combined with dynamical models, 
the data suggest that the zodiacal dust is largely cometary, 
rather than asteroidal, in origin. 



1. INTRODUCTION 



The motion of interplanetary dust particles contains im- 
portant information about their origin, distribution, and 
orbital evolution. At optical wavelengths, dust with radii 
« 10-100/im that lie within « 3 AU of the Sun scatters 
the incident solar radiation to produce zodiacal light, and 
the relative motion of the dust modifies the location and 
shape of solar spectral lines [1-3]. The fraction of zodi- 
acal dust with a cometary or asteroidal origin is not well 
constrained at present [4-5], and the kinematics of these 
two components may shift the velocity and widths of the 
spectral features in unique ways. However, the low sur- 
face brightness of zodiacal light has, until recently, lim- 
ited the observability of this effect, requiring a combina- 
tion of high sensitivity and high spectral resolution. Pre- 
vious work has shown the dust to be on prograde orbits, 
but different groups reported uncertain and contradictory 
results regarding the details of the orbital properties of 
the dust [6-7]. Here, we report on new measurements 
of scattered solar Mg I A5184 absorption line in zodiacal 
light with the Wisconsin H- Alpha Mapper (WHAM), and 
compare the observations with predictions from dynami- 
cal models of the zodiacal dust cloud. 



WHAM consists of a 15cm, dual-etalon Fabry-Perot 
spectrograph coupled to a 0.6m siderostat, and produces 
an average spectrum over a 1 ° circular field of view with 
a 12 km/s resolution within a 200 km/s spectral window. 
It is located at the Kitt Peak National Observatory in Ari- 
zona, and is entirely remotely operated. It is specifically 
designed to study extremely faint, diffuse optical light at 
high spectral resolution [8]. We have recorded spectra 
centered on the Mg I line at 5 183.6A toward 49 directions 
along the ecliptic equator, with two directions at high 
ecliptic latitude on the nights of 2002 November 4 and 
5 [9]. The unprecedented capabilities of WHAM allowed 
us to identify and remove several weak atmospheric emis- 
sion lines that probably affected the results of previous 
investigations [6-7]. The terrestrial lines are stronger on 
the red side of the Mg I line. If they are not accounted 
for, they shift the line centroid to negative velocities and 
produce asymmetries in the line profile shape. 

A sample spectrum demonstrating the high resolution 
and sensitivity of the observations is shown in Figure 1 . 
This spectrum was taken toward the north ecliptic pole, 
with the brightness of the line given in units of milli- 
Rayleighs per km/s. The shaded solid line is a spec- 
trum taken at twilight, where the light is scattered off the 
Earth's atmosphere which is at rest relative to the zodia- 
cal cloud. The two weak absorption lines in the twilight 
spectrum near H-50 km/s are other solar lines of Fe I and 
Cr I. The line centroid, width, and area of each spectrum 
were measured using a least-squares fitting technique, 
with the twilight spectrum used as a template. Figure 2 
shows the change in velocity centroid with solar elonga- 
tion (helio-ecliptic longitude, e) for directions along the 
ecliptic equator. The data show a clear prograde kine- 
matic signature. We see no evidence for a net radial out- 
flow (at e=180°) or an east/west asymmetry as reported 
by previous authors [6-7]. 

Figure 3 shows the change in the full-width at half- 
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Figure 1. Spectrum of the zodiacal light taken toward 
the north ecliptic pole, centered on the Mg I A5184 ab- 
sorption line. The solid line shows the Mg I line from 
the twilight sky, and represents an unperturbed reference 
spectrum. Note the flat-bottomed profile and enhanced 
line width relative to the twilight spectrum. This suggests 
that a significant population of the zodiacal dust particles 
are on non-circular orbits. 
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Figure 2. Velocity centroid of the Mg I absorption line as 
a function of helio-ecliptic longitude. The observational 
data are shown as solid circles, with model predictions 
are shown as lines. The data indicate the dust follows 
prograde orbits, with no significant net radial outflow or 
east/west asymmetry. The observations are better fit by 
models in which the dust follows elliptical or cometary- 
like orbits (see text). 



3. COMPARISON TO MODELS 



maximum of the lines as a function of elongation for di- 
rections along the ecliptic equator. The solid horizontal 
line near 58 km/s is the width of the unperturbed twilight 
spectrum. The Figure shows that along the ecliptic plane, 
the line is broadened by 10-20 km/s relative to the so- 
lar line, with no significant trends with elongation. The 
large intrinsic width of the solar line makes these mea- 
surements difficult to quantify accurately and contributes 
to the large uncertainty compared to the velocity centroid 
data. 

If all the dust particles were on pure circular orbits cen- 
tered on the Sun, the zodiacal and twilight line profiles 
would be nearly identical toward the anti-solar direction. 
The broadening of the zodiacal profiles can only be ex- 
plained by particles having radial components to their or- 
bital motion. The data imply that a significant fraction of 
the dust follows elliptical orbits. 

Furthermore, the substantial broadening of the profile to- 
ward the ecliptic pole (Figure 1), which is even greater 
than that toward the anti-solar direction (see [8]), implies 
a population of particles that have significant components 
of their orbital velocities projected perpendicular to the 
ecliptic plane. For an orbital speed of 30 km/s near 1 
AU, and neglecting the radial motion associated with ec- 
centricity, the approximately ±15-20 km/s broadening at 
the base of the profile suggests a distribution of inclina- 
tions extending up to ~30°-40° with respect to the eclip- 
tic plane [9]. 



The relationship between the observational data and the 
kinematic properties of the zodiacal dust cloud is a clas- 
sical inversion problem. The shape of the observed line 
profiles is determined by the population of dust particles 
of varying size, radial distance, scattering function, and 
relative motion along the line of sight. To infer the orbital 
properties of the particles that comprise the dust cloud, 
the data need to be compared to a range of predictive dy- 
namical models. There has been considerable work done 
in the field of modeling the zodiacal dust cloud that take 
into account a wide range of observational characteristics 
[e.g., 10-12]. However, few models explicitly consider 
the dynamics of the cloud that can be directly compared 
to our data. Below we discuss a few models from the 
literature, including some new models, that estimate the 
change in both the velocity centroid and line width with 
elongation. 

Some of the models that compare favorably to the data 
are shown in Figures 2 and 3. The solid black line is a fit 
to a model from [13]. This model describes particles on 
prograde, elliptical orbits with eccentricities uniformly 
distributed between and 1, with randomly distributed 
perihelions. Their model did not include the influence of 
radiation pressure, and the particles were confined to the 
ecliptic plane. The model fits the centroid data well (Fig- 
ure 2), but strongly overestimates the width of the lines 
(Figure 3). The inclusion of radiation pressure and/or in- 
clined orbits could provide a better match to the observa- 
tions [9,14]. 
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Figure 3. Full-width at half-maximum of the Mg I ab- 
sorption line as a function of helio-ecliptic longitude. The 
observations are shown as solid circles with the model 
data shown as various lines. The solid horizontal line rep- 
resents the width of the absorption line taken at twilight. 
The motion of the zodiacal dust broadens the line by « 
20 km/s relative to the twilight spectrum, with no sig- 
nificant changes in line width with elongation. With the 
exception of the model of [13], the models are generally 
well-matched to the observations, within the uncertainty. 



The colored lines in Figure 2 and 3 are models from [15- 
16], which trace the motion of different populations of 
dust particles subject to gravity, radiation pressure, and 
drag forces. The individual dashed and dotted lines rep- 
resent particles with asteroidal and various cometary tra- 
jectories, with a ratio of radiation pressure to gravitational 
force of /? = 0.002. Some models considered particles 
with different /3-values, and some with trans-Neptunian 
orbits, but those models are generally poorly matched to 
the data and are omitted for clarity. In Figure 2, we see 
that a good match is provided by the cometary particles 
on inclined, eccentric orbits compared to the asteroidal 
particles. In Figure 3, all of the models fall within the 
large scatter in the data, do not allow us to discriminate 
between the different models. Future observations of in- 
trinsically narrower lines, such as Fe I, will aid in using 
line widths to assess the asteroidal and/or cometary-type 
orbits of the zodiacal dust. 



4. SUMMARY AND FUTURE WORK 



Observations of scattered solar absorption lines in the zo- 
diacal light are a powerful technique for exploring the 
kinematics of the zodiacal dust cloud. Our data are fit 
well by models that contain particles on elliptical orbits 
that are inclined to the ecliptic plane. This suggests that 
most of the dust in the zodiacal cloud has a cometary ori- 
gin, in agreement with [4-5] . The interpretation of our ob- 
servations is highly model dependent, and we emphasize 



the importance of exploring a range of predictive mod- 
els in assessing our conclusions. Higher signal-to-noise 
observations covering a larger fraction of the ecliptic sky 
that include other, more intrinsically narrow, absorption 
lines will provide a more complete picture of the kine- 
matics of the zodiacal dust cloud. New dynamical mod- 
els that investigate a wider range of dust properties and 
orbital parameters can provide strong, quantifiable con- 
straints on the nature of the zodiacal dust cloud when 
compared to these new observations. 
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